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INTRODUCTION 

Analysis of interaction occurring between space de- 
bris and orbiting structures is of great interest to the 
planning and survivability of space assets. Computer 
simulation of the impact events using hydrodynamic 
codes can provide some understanding of the pro- 
cesses but the problems involved with this fundamen- 
tal approach are formidable. First, any realistic simu- 
lation is necessarily three-dimensional, e.g.. the 
impact and breakup of a satellite. Second, the thick- 
ness of important components such as satellite skins 
or bumper shields are small with respect to the dimen- 
sion of the structure as a whole, presenting severe zon- 
ing problems for codes. Thirdly, the debris cloud pro- 
duced by the primary impact will yield many secondary 
impacts which will contribute to the damage and pos- 
sible breakup of the structure. Characterization of the 
debris cloud requires accurate fragmentation model- 
ing as well as accurate tracking of the fragments 
through large regions of void. For these reasons hy- 
drodynamic simulation of hypervelocity impact and 
breakup of orbiting structures is extremely difficult. 
We have approached the problem by choosing a rela- 
tively new computational technique that has virtues 
peculiar to space impacts. The method is called 
Smoothed Particle Flydrodynamics (SPFI). In this pa- 
per we describe the SPFI method and why we believe 
that it can be used to answer many questions concern- 
ing the survivability of space assets due to kinetic im- 


pacts. We also present several calculations to show the 
power of SPH towards such problems. 

SPH BASICS 

Smooth Particle Hydrodynamics is unique in 
computational fluid dynamics in that SPH uses no 
grid. It was the genius of the inventors, Lucy (1977) 
and Gingotd & Monaghan (1977) to figure out how to 
take a derivative (get the force on a fluid element) 
without using a mesh. Previously a mesh was the only 
known way to compute a spatial derivative using fi- 
nite differences. The mathematical theory of SPH 
will not be discussed here. The reader is referred to 
Gingold & Monaghan (1977,1982), Monaghan 
(1982,1985) and Monaghan & Gingold (1983) for de- 
tailed treatment of the subject. We present here only 
some basic features as discussed by Benz (1989) that 
are necessary to understand the method. Consider a 
function/, a kernel W which has a wadth measured by 
the parameter h, and the following equation: 

<flr) >= \ w{T-x\hwW. 0) 

If w r e impose a normalization condition such that the 
integral of W is unity, then it follows that 

< Ar) > h ~*° A T ) ■ (2) 
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Relation (1) therefore defines the kernel estimate <f> 
of /. If W is the Dirac delta function then we have the 
equality </>=/. Now suppose / is known only at N 
discrete points that are spatially distributed according 
to the number density distribution: 

N 

n(r) = ^<5(r-r y ) (3) 

;= i 

If the number density at Tj is written as 


< n(r,) > = 


P(*>) 


(4) 


thus introducing the concept of particle mass (m), the 
following equation can be derived: 


<f(r) > =J j f J W(r-r J .h)^- (5) 
j 6 > 

This equation defines a procedure for transforming 
integral equations to particle equations and is 
therefore called "integral evaluation by the particle 
method." The choice of kernel or "smoothing 
function" is discussed by Monaghan and Lattanzio 
(1985). The Wmost frequently used in SPH codes is 
a B-spline with compact support which goes to a zero 
at a distance 2h from its peak. When the conservation 
laws of fluid dynamics are cast into the SPH 
framework using the procedure described above, the 
following equations are obtained. 

q, = ^ m,W,j ( 6 ) 
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Equation (6) is the density computation of particle T' 
using the masses of neighboring "j" particles. The ac- 
celeration of particle "P is given by (7) and the evolu- 


tion of the specific internal energy (e) is described by 
equation (8). These equations also involve the pres- 
sure (P) and the artificial viscous pressure n. Terms 
involving material strength are omitted here but are 
discussed by Libersky and Petschek (1991). It should 
be emphasized that equations (6) thru (8) are the con- 
tinuum equations of fluid dynamics cast into a discrete 
Lagrangian frame by kernel interpolation. Therefore, 
the Smoothed Particle Hydrodynamics method is a 
true and complete hydrodynamic calculational proce- 
dure. 

IMPORTANCE TO SPACE IMPACTS 

Simulating the mechanics of irreversible processes 
that take place between a structure and a projectile 
during a hypervelocity collision in space is extremely 
difficult. There are three main reasons for this. First, 
problems are three-dimensional. Second, the thick- 
ness of important components such as satellite skins, 
electronic components or bumper shields are small 
with respect to the dimension of the structure as a 
whole. Three-dimensional calculations using mesh- 
based codes do not seem feasible for such problems. 
Thirdly, the debris cloud produced by the primary im- 
pact produces secondary and tertiary impacts impor- 
tant to the overall damage and breakup. Simulation 
of these events requires detailed characterization of 
the debris cloud which in turn requires good fragmen- 
tation modeling as well as accurate tracking of the 
fragments through large regions of void. Eulerian 
codes have difficulty tracking sub-grid scale frag- 
ments through the mesh. Also, large regions of void 
within the structure need to be zoned in anticipation 
of material arriving there at some later time. SPH suf- 
fers from neither of these difficulties because there is 
no mesh. Following the debris cloud through large 
regions of void presents no difficulty to SPH. Further- 
more, interfaces between materials in a problem con- 
sisting of several materials are accurately tracked. 
These virtues of SPH are true of any Lagrangian code, 
except that mesh-based Lagrangian codes cannot 
treat large fluid distortions. Obviously, hypervelocity 
impact produces highly distorted flows. It can be said 
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that SPH contains the best features of grid-based Eul- 
erian and Eagrangian methods without the limitations 
of either. The price to be paid is that SPH is slower 
than the other methods, (having to determine new 
neighbors each computational cycle) but not prohibi- 
tively so. The efficiency of SPH appears to be much 
better than Eulerian in 3D (Durisen, 1988) and we are 
paying close attention to vectorization and paralleliza- 
tion. The ability of SPH to accurately track debris re- 
sulting from hypervelocity impact through large re- 
gions of void in 3 dimensions and compute the impact 
of the debris on other parts of the structure make the 
method extremely attractive for space applications. In 
working towards that goal we have performed several 
two-dimensional calculations to evaluate SPH and get 
a feel for how impact events might damage and lead 
to breakup of the large-scale structure. 

SYNTHESIS OF A PREDICTIVE BREAKUP 
MODEL 

Determination of processes that contribute to the fail- 
ure and ultimately to the breakup of a complex and in- 
tegrated space asset under impulsive loading is of vital 
interest. Currently, such a predictive model does not 
exist. Until recently, virtually all breakup modeling of 
a spacecraft under intensive loading used a phenome- 
nological approach. Although this approach provides 
some qualitative measure of the interactions, it lacks 
the physics necessary to identify those processes that 
control generation of the resulting debris cloud. To 
characterize the environment of a debris cloud accu- 
rately (in terms of debris mass, velocity and number 
distribution), a first principles physics based predictive 
model has been synthesized. This approach considers 
the synthesis of a total predictive model based on the 
response of elementary components. Internal compo- 
nents of a generic satellite are shown in Figure 1 and 
an example matrix of calculations leading to the total 
predictive model is shown in Figure 2. The design of 
these and other calculations is motivated by the vari- 
ous geometries of the components internal to the sat- 
ellite. Two simulations of local impact events on small 
regions of a large scale structure are shown in Figure 


3 and Figure 4. These calculations give us detailed un- 
derstanding of how the debris cloud interacts with 
nearby structural elements producing damage. 

Figure 3a is a particle plot showing initial conditions 
for the impact of an tantalum projectile on an alumi- 
num (2024-T86) frame. Each arm of the frame had 
length 20.5 cm and thickness 0.5 cm. The projectile 
had 1.0 cm sides, a speed of 7 km/s and a 60 degree 
impact obliquity. The calculation was performed in a 
two-dimensional Cartesian frame of reference. A 
Gruneisen equation of state and an elastic-perfectly 
plastic constitutive model were used to describe the 
metals. The calculation used 12,240 particles and the 
smoothing length was 0.08333 cm. Results of the cal- 
culation at 1 and 20 microseconds are shown in Fig- 
ures 3b and 3c respectively. A very large opening is 
created in the first plate impacted due to the large im- 
pact obliquity. This "hole" is approximately 7 times 
the initial projectile size. One end of the plate is bent 
inward. There is a large splash of material moving up- 
ward and away from the structure typical of high speed 
cratering events and a debris cloud expanding towards 
the other arm of the structure. This cloud has frag- 
ments of various sizes but most of the mass is concen- 
trated in the part of the cloud that is about to impact 
the second structural arm. 

Results of a similar calculation employing a slightly 
more complex aluminum structure are displayed in 
Figure 4 where particle plots at four different times (0, 
10, 20, 30 jis)are shown. In this calculation the debris 
cloud, resulting from the impact upon one member of 
the aluminum frame, is seen impacting other structur- 
al elements. Severe damage is seen on the upper pan- 
el. In fact, a secondary debris cloud has been pro- 
duced by the interaction. It is easy to see from this 
calculation, that in an actual satellite with many more 
composite parts, how a cascade of debris clouds could 
form to cause massive breakup of the structure. Some 
damage is also seen on the outside panel furthest from 
the impact. These calculations were performed on a 
1 megaflop machine and required approximately 1 
hour of cpu time per 10 ps of simulation time. 
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Target Name: chucksat.tp4 
Aimpoint(X»Y,Z) (cm) = ( 0.00 0.00 45.50) 

Observer Azimuth (degrees) = 91 .00 

Observer Elevation (degrees) = 0.00 



Figure 1. A Generic Satellite 
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Figure 1A. Interior Components of Satellite 
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Matrix of Calculations Leading to a 
Satellite Breakup Model 
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Figure 2, Matrix of Calculations 
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il configuration. 


after initial impact 


lum projectile penetrating an aiitr.i 
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Material Plot Material Plot 
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Figure 5. Smoothed Particle Hydrodynamics simulation of a liquid filled 
steel pellet impacting multiple spaced plates at 2 km/s and high obliquity. 
The first plate impacted is aluminum amd the remaining two are steel. The 
calculation was done in two-dimensional Cartesian coordinates using the 
SPH code MAGI. Frames a - d show results at 20, 180, 325, and 647 ps. 
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Figure 6. Smoothed Particle Hydrodynamics simulation of a liquid filled 
aluminum tank hit by an aluminum projectile moving at 7 km/s. The 
calculation was done in two-dimensional Cartesian coordinates using the 
SPH code MAGI. The liquid particles are gray-scaled by density and the 
aluminum particles are shown in black. Frames a - d show results at 3, 30, 
60, and 120 ps. 
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Several calculations have been presented - two hyper- 
velocity impacts on simple structures, and two other 
impact problems possessing interesting features. We 
are encouraged by the results. The code has been ex- 
tended to three-dimensions and we are currently per- 
forming test problems. 
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ADDITIONAL IMPACT CALCULATIONS 

Two additional calculations are presented in order to 
show' the power of SPH towards impacts. Figures 
5a-5d are particle plots showing a liquid filled projec- 
tile impacting a series of spaced plates at high obliqui- 
ty. The projectile casing w'as steel and the simulated 
liquid was modeled as carbon-tetrachloride. The 
length of the projectile was 3.75 cm and the thickness 
2.50 cm. The first plate to be impacted was aluminum 
with thickness 0.254 cm, and the remaining two plates 
were steel with thicknesses 0.48 cm and 1.58 cm re- 
spectively. The projectile speed was 1.6 km/s. The 
problem was run in a two-dimensional rectangular 
Cartesian frame of reference (plane-strain) with ap- 
proximately 4,000 particles. The smoothing length was 
0.127 cm. An elasctic-perfectly plastic constitutive 
model and Gruneisen equation of state were used to 
describe the solids. No explicit failure model was in- 
cluded in the calculation. Figure 5a shows the alumin- 
ium plate being impacted (t = 20 gs). Figure 5b shows 
the calculational result at 180 ps. I he projectile has 
struck the first steel plate and the low^er aluminum 
plate which was set in motion by the impact has also 
hit the steel and started to buckle. At 325 ps, as shown 
in Figure 5c, the aluminum plate has continued to 
buckle and the portion of the broken steel plate at the 
bottom of the figure has rotated as a result of the im- 
pact torque. These are two interesting features of the 
strength model operating in the code. It is very en- 
couraging to see these effects captured by the code. 
Other interesting features are the crack formation at 
the back of the thick steel plate and the enhanced up- 
ward momentum of the first steel plate due to the cra- 
ter splash from the projectile on the second plate. We 
expect these kinds of secondary effects to be impor- 
tant in the actual satellite breakup from impact. An 
obvious key feature of simulation is the "plug" of steel 
plate seen in Figure 5d. This is not due to any explicit 
failure model in the code but results from intrinsic 
model fracture in SPH. In response to elastic waves 
in the plate, some particles find themselves outside of 
the smoothing length range of communication with its 


neighbors. This separation manifests itself as a 
"crack" which propagates through the metal. We can- 
not, at this point, claim that this effect is actually a 
physical one. Nevertheless, it appears that the code is 
trying, on its own, to accommodate failure, and the re- 
sults seem physical. This calculation required 2 
hours of cpu time on the CRAY 2 at the Phillips labo- 
ratory. 

Results of a hypervelocity impact (7 km/s) simulation 
of an aluminum projectile into a water-filled alumi- 
num tank are shown in Figure 6 where the particles are 
gray-scaled according to their density. The diameter 
of the tank was 25 cm initially and the wall thickness 
was 0.5 cm. The projectile length was 3.0 cm and the 
thickness was 2.0 cm. The geometry was 2D Carte- 
sian. Approximately 20,000 particles were used in the 
simulation. The smoothing length was 0.2 cm. Notice 
the shock propagation in the water and along the tank 
inner surface which comes to a focus at the rear of the 
cylinder and then strongly reflects back. The reflected 
shock is evidenced by the flattening of the particle dis- 
tribution just ahead of the projectile. Notice also the 
rapid deceleration of the aluminum projectile and its 
large deformation. I his calculation took 7 hours on 
the CRAY 2. For these two specialized problems 
there is no experiment data to which the simulations 
can be compared, so we must be careful not to draw 
unjustified conclusions about the codes performance. 
However, we can get a feel tor how the code responds 
to difficult impact problems with the goal of extending 
the calulations to three-dimensions with extensive 
comparison to experiments. 

CONCLUSIONS 

The Smoothed Particle Hydrodynamics (SPH) ap- 
proach to computational fluid dynamics has been 
briefly described with emphasis on its natural ability 
to model hypervelocity impact on orbiting structures. 
Our goal is to exploit these virtues of SPH towards the 
development of a complete structural breakup model 
in order to answer important questions concerning the 
survivability of space assets due to kinetic impacts. 
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(4a) 


Initial configuration. 


Figure 4. 



(4b) 10 microseconds after initial impact. 


Material plot of a tantalum projectile penetrating an aluminum 
A-frame structure at 60° obliquity. 
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